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1. Introduction 

W. Bietenholz has explained our procedure for extracting Wilson coefficients from lattice mea- 
surements in his proceedings |JXJJ . In this report we will show how these coefficients can be used 
to reconstruct Compton scattering amplitudes, and to extract their Nachtmann moments. We start 
with a toy example, looking at the Nachtmann moments of off-shell quarks in Landau-gauge back- 
ground fields. 

The calculations we report here are done with overlap fermions on a quenched background. 
Overlap fermions were chosen for their superior chiral symmetry properties. The negative mass 
parameter is p = 1.4. We used a 24 3 x 48 lattice, with a lattice spacing a = 0.095 fm. The results 
shown here all have bare mass am q = 0.028. We have always taken the scattering momentum q 
along a lattice diagonal, qoc (1,1,1,1) for maximum lattice symmetry. We employ three different 
values for the magnitude of a 2 q 2 , so that we can begin to investigate scaling in q 2 . All the Green's 
functions have been 0(a) improved using the prescription of The electromagnetic current 
Jjx(x) is represented by the local current \jf(x)Yu.Y(x), with overlap 0(a) improvement. 

2. Operator Product Expansion 

We express the electromagnetic scattering tensor W^ v as a sum over local operators (the Oper- 
ator Product Expansion or OPE). In each term we have a separation of scales, all dependence on the 
quark momentum is in the matrix element (\\r(p)\J ^(q)jl(q)\\\r(p)) , all dependence on the photon 
scale q is in the Wilson coefficient C™ v (q). At present we are only considering flavour non-singlet 
processes, so we do not include any purely gluonic operators in the sum, 

w, v ( P ,q) = (w(p)\JMA{q)W(p)) = ^c;M{w(p)W m W(p)) ■ CD 

m 

We include quark bilinear operators G m with up to 3 covariant derivatives in this sum. When all 
possible Dirac structures are taken account of, there are potentially 1360 different operators, and 
1360 Wilson coefficients C™ v , in the sum. We reduce this number by exploiting lattice symme- 
tries. We choose q along a lattice diagonal, i.e. q oc (1,1,1,1). With this choice there are only 67 
independent Wilson coefficients in the expansion of a diagonal element of W, such as W44 ph. 

W^ v (p,q) is a fairly complicated object, it depends on p,q,/x,v and on the Dirac indices of 
the incoming and outgoing quark. 

We simplify by just looking at unpolarised quarks (later, we plan to analyse spin-dependent 
quantities too). Taking the trace 

Tnv(p,q) = ^Tr{S- 1 (p)W flv (p,q)} (p 2 +m 2 ) , (2.2) 

where 5 _1 is the inverse quark propagator, removes the Dirac-index structure. The propagator 
cancels all Z v factors, the only renormalisation we need is a factor of Z 2 , to correct for using local 
currents for . 

If we consider unpolarised quarks there are 4 tensor structures which can occur in the scattering 
tensor: 

T^v = Sf! V Wi + p^iP v W2 + (p^q v +q ti p v )W 4 + q^q v W 5 . (2.3) 
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(W3 and W(, are reserved for structures possible in neutrino scattering and the polarised target case.) 
The Wj form factors can only depend on invariants q 2 , p ■ q, p 2 . When we consider scattering on 
physical hadrons, we can use electromagnetic gauge invariance to reduce the expansion from four 
terms to two, namely F\ and F2. When however we consider an off-shell quark this argument no 
longer applies, and all four structures are independent. 

As a first simple case we consider the polarisation trace of T^ v , 

Tfijj, = 4Wi + p 2 W 2 + 2p ■ qW 4 + q 2 W 5 . (2.4) 

Advantages of this choice are that averaging over the direction of (the photon polarisation) 
simplifies the rotation group theory considerably, and that this quantity only involves diagonal 
elements of W, which we have analysed more completely. (We have gathered data on off-diagonal 
components, /I 7^ v, but this has not yet been fully analysed.) Summing over all polarisations 
also increases the symmetry, there are only 22 independent Wilson coefficients in the OPE of 7^, 
which is a considerable reduction compared to the 67 coefficients needed for a single diagonal 
component such as T44. 

In [Q] Nachtmann proposed some quantities (the Nachtmann moments) with particularly sim- 
ple Operator Product Expansions. Nachtmann considered scattering from on-shell targets, but we 
are interested in off-shell targets, so we have to generalise the formulae in [0]. The Nachtmann 
moments, /!„, are defined by splitting T up into components of definite spin, n, 

T^(p-q,q 2 ,p 2 ) = 2 £ (4) ' ^( cos ^M^V) ( 2 - 5 ) 

n even V 9 / 

where 6 is the angle between p and q, 

p ■ q = \p\ \q\ cosd (2.6) 

and U n is a Chebyshev polynomial of the second type, [^J. The U n are the 4-dimensional equivalent 
of the familiar 3-dimensional spherical harmonics. We can use orthogonality of the Chebyshev 



polynomials to project out single Nachtmann moments from (2.5), 



J o — sin 2 6> C/„(cos0) T^(p-q,q 2 ,p 2 ) = IX n (q 2 ,p 2 ) . (2.7) 

The Euclidean integral ( |2"77| ) involves real values of p • q in the range — \p\ \q\ < p q < \p\ \q\. 

Following [Q] we can relate the integral to Minkowski physics by allowing p ■ q to become 
complex, while keeping both p 2 and q 2 real and positive, see Fig.|]. Ignoring the possible com- 
plications due to confinement, we expect the amplitude T to have branch-points at the thresholds 
for producing on-shell quarks, at p ■ q = ±^(q 2 + p 2 + m 2 ), with cuts reaching out to infinity. We 
define the discontinuity across these cuts by 

2niD(p-q,q 2 ,p 2 ,m) = T^(p-q + ie,q 2 ,p 2 ) -T^(p-q-ie,q 2 ,p 2 ) . (2.8) 

For colour singlet hadrons this discontinuity is a physically measurable total cross-section, but 
of course the cross-section for deep inelastic scattering on an off-shell quark target is something 
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Figure 1: 77ie complex p ■ q plane, for p 2 ,q 2 both positive (space-like). The Euclidean integral (2.7) for the 
Nachtmann moment \i n runs between p ■ q = — \p\ \q\ and p ■ q = +\p\ \q\, the Minkowski integral (2.9) runs 
along the cut from the threshold at p -q = j(q 2 + p 2 + m 2 ) to °°. 



we can only measure as a Gedankenexperiment. Assuming that T is an analytic function of the 
complex variable p-qwe, can write down dispersion relations which give the result of the Euclidean 
integral (2.7) as an integral involving the discontinuity (valid for even n,n> 2). 

D{p-q,q\p 2 )(q 2 f 



Vn(q ,P )=2 j_ d{p-q) 

' [p-q+V(p-q) 2 -p 2 q : 

where & is the threshold for the production of on-shell particles, 

& = Uq 2 +p 2 + m 2 ) . 



n+l 



(2.9) 



(2.10) 



In general the Minkowski integral (2^9) is complicated, but in the Bjorken limit, when {p ■ q) S> 
p 2 q 2 , we can simplify (2.9) by changing to the integration variable 

-.2 



q 



2p-q 



(2.11) 



Eq. ( p.9| ) becomes 

Vn(q 2 ,p 2 ) 



dx 



x n - x D{p-q,q 2 J p 2 ) 

n+l 



+ ■ 



dxx* 1 1 D(p ■ q,q l ,p z 



(2.12) 



and we see that the Nachtmann moment tends to a simple x moment. 

In the rest of this report we concentrate on the smallest interesting spin, n = 2. From ( 2.12 ) 
we see that at large q 2 , jl2 corresponds to (jc). 

In our discussion of Nachtmann moments we have assumed full rotation symmetry. However, 
lattice operators are classified under the hypercubic group, and will normally be mixtures of repre- 
sentations of the full Euclidean rotation group. The 22 operators present in the expansion of 
fall into three hypercubic classes 
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spin 0, spin 4, + higher 8 operators 

spin 2, spin 4, + higher 13 operators 

spin 4, + higher 1 operator 

For our first look at spin 2, we simply keep all the operators in the middle group, and discard the 
others. The leading spin 2 operator (for qoc (1,1,1,1) ) has the form 

Y.WiynDy + YvD^w (2.13) 

i.e. it is a symmetric, off-diagonal tensor. 

We now reconstruct a projected spin 2 amplitude, [r^] jp! -„ 2> by taking the product of the 
matrix elements of the 13 spin 2 operators with their Wilson coefficients, as determined in [jj]]. 



Figure 2: The 4-dimensional spherical harmonic p 2 q 2 U2(cos 9) — 4(p ■ q) 2 — p 2 q 2 - The main lobes point 
along the directions ±q. There are nodes when p is 60° and 120° from q; the equatorial 'donut' has the 
opposite sign to the main lobes. 

Because we have averaged over quark spin and photon polarisations, the only direction left in 
our problem isgoc (1,1,1,1) and the only spin-2 4-d spherical harmonic that contributes is 

p 2 q 2 U 2 (cosd) = 4(p-q) 2 -p 2 q 2 , (2.14) 

which is illustrated in Fig.||. When we put in the fact that q oc (1,1,1,1) 

4{p-q) 2 -p 2 q 2 ^2(p iP 2 + piP3 + ---+P3P4)q 2 ■ (2.15) 

Spin 2 operators, such as ( |2.13| ), should have an expectation value proportional to (2.15). 



We can find the \i 2 for the quark from (2.7), 




(2.16) 

so if we plot q 2 [T^^\ spjn2 against 2 (L jU<v P\iPv) we should see a straight line passing through the 
origin, with a slope equal to \L 2 - 

We have data for three different q values, and from 15 to 32 p values (depending on q), chosen 
to give a good coverage of directions, so we can see whether the projected amplitude really follows 
the spherical harmonic, and find jj.2 for a quark. 
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3. Results 




-0.2 0.2 0.4 0.6 0.8 



Figure 3: The trace q 2 \T^fi\ spin2 plotted against 2 (^ <v P/iPv), in lattice units. These points should fall 
on a straight line through the origin, with a slope proportional to the Nachtmann moment corresponding 
to (x), see eq.(2.16). The black squares are data from q 2 — n 2 / (2.25a 2 ) ms 19 GeV 2 , blue circles, q 2 — 
n 2 /(4a 2 ) w 10.6 GeV 2 , red triangles q 2 = n 2 /(9a 2 ) « 4.7 GeV 2 . 

In Fig.|| we show the spin-projected Compton amplitude plotted against 2 (L^< v P^Pv)- The 
quantity 2 (Y,n<v PuPv) is positive for momenta in the main lobe of the harmonic in FigfJ, zero 
for momenta in the node, and negative for momenta in the equatorial ring. 

At each q 2 value the amplitudes follow a straight line rather closely. Several quarks have 
momenta at exactly 60° angle to q, these fall on the node of the spherical harmonic, (zero on the 
horizontal axis), and they have amplitudes very close to 0, as predicted. Data from all three q 
values are plotted, we see that they scale fairly well, q 2 changes by a factor of 4, from 4.7 GeV 2 to 
19 GeV 2 . 

Rough values for the Nachtmann moment \i2 (which is a measure of (jc) at a scale ~ q 2 ) are 

q 2 = 4.7 GeV 2 \l 2 = 0.44(9) 

q 2 = 10.6 GeV 2 ^2=0.73(5) 
q 2 = 19 GeV 2 \i 2 = 0.60(5) 

These values probably include fairly large lattice artefacts oc a 2 q 2 . We will try to correct for these 
by looking at tree-level lattice artefacts, which might lead to significant numerical changes. 
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4. Conclusions and Prospects 

We have seen how we can project out spin components for Nachtmann moment operator ex- 
pansions. In the channel we looked at (spin 2), simply filtering on the lattice symmetry seems to 
work fairly well, we don't see any sign of spin 4 contamination distorting the straight-line behaviour 
ofFigf 

We don't see any strong dependence of /I2 on the quark virtuality p 2 , this may be a little 
unexpected. 

There are several more things we could do. Here we have concentrated on the Nachtmann 
moment coming from the polarisation trace W^, we should also look at the moments for the other 
components of W^ v . Particularly interesting in the context of off-shell quarks would be to look 
at the OPE for q^W^qv. This should give a combination of operators which can be non-zero for 
off-shell quarks, but which should vanish on-shell. The three-point functions for this combination 
of operators should show contact terms, but no plateau or exponentially decaying terms. 

The antisymmetric parts of the scattering tensor, W^ v — W Vjl contain the information needed 
to investigate spin-dependent structure functions. We have collected the data needed for this calcu- 
lation. 

From looking at our problem in tree level we see that 0(a 2 q 2 ) artifacts may be important in 
some channels. We want to investigate these artifacts further, and see whether we can use tree-level 
results to reduce their severity. 
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